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We examine simulations of core-collapse supernovae in spherical symmetry. Our model is 
based on general relativistic radiation hydrodynamics with three-flavor Boltzmann neutrino 
transport. We discuss the different supernova phases, including the long-term evolution 
up to 20 seconds after the onset of explosion during which the neutrino fluxes and mean 
energies decrease continuously. In addition, the spectra of all flavors become increasingly 
similar, indicating the change from charged- to neutral-current dominance. Furthermore, 
it has been shown recently by several groups independently, based on sophisticated su- 
pernova models, that collective neutrino flavor oscillations are suppressed during the early 
mass-accretion dominated post-bounce evolution. Here we focus on the possibility of col- 
lective flavor flips between electron and non-electron flavors during the later, on the order 
of seconds, evolution after the onset of an explosion with possible application for the nu- 
cleosynthesis of heavy elements. 

1 Introduction 

Explosions of massive stars are related to the formation of a shock wave, which forms when 
the collapsing stellar core bounces back at nuclear matter density. During collapse, the stellar 
core dclcptonizes so that a low central proton-to-baryon ratio, given by the electron fraction of 
Y e ~ 0.3, is reached at bounce. The conditions obtained at bounce depend sensitively on the 
weak interaction scheme and the equation of state used. Fig. [1] illustrates the radial evolution 
of selected mass elements. Before bounce, the infalling mass elements correspond to the central 
iron-core while the outer layers of the progenitor are basically unaffected from the central 
happenings. After bounce, the shock wave propagates outwards and stalls on a timescale of 
5-20 ms due to energy losses from heavy-nuclei dissociation and ^-escapes emitted via large 
numbers of electron captures during the shock passaged across the neutrinosphercs. As a result 
of energy loss, the expanding dynamic bounce shock turns into a standing accretion shock 
(SAS). For the early shock propagation and the position of the zy e -sphere, see the red solid and 
magenta dash-dotted lines in Fig. [1] The post-bounce evolution is given by mass accretion 
onto the SAS and neutrino heating, dominantly via v e and v e absorption at the dissociated free 
nucleons, behind the SAS on timescales on the order of 100 ms. 

Several explosion mechanisms have been explored; the magneto- rotational [1] , the dumping 
of acoustic energy [5] and the standard scenario due to neutrino heating [3] . Recently, it has been 
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Figure 1: Sketching the 
evolution of selected mass 
elements during core col- 
lapse, bounce, post-bounce 
accretion and onset of ex- 
plosion. Color coding 
is according to the domi- 
nant composition (light red: 
heavy nuclei, blue: 4 He, yel- 
low: light nuclei and free 
nucleons). The solid red 
and dash-dotted magenta 
lines mark the positions of 
shock and ncutrinosphere. 
The dashed lines mark the 
evolution of interfaces be- 
tween different composition 
layers of the progenitor. 

shown that a quark-hadron phase transition can lead to the formation of an additional shock 
wave that can trigger explosions [4] [5]. In this article, we explore standard neutrino-driven 
explosions in spherical symmetry of the low-mass 8.8 M Q O-Ne-Mg-core and more massive 
iron-core progenitors. For the latter, where neutrino-driven explosions cannot be obtained 
in spherical symmetry, we enhance neutrino heating in order to trigger explosions. Fig. Q] 
illustrates the standard neutrino-driven explosion of a 15 M Q progenitor, for which the accretion 
phase ends at about 450 ms post bounce with the onset of explosion. The SAS turns into 
a dynamic shock which expands continuously to increasingly larger radii (see the red solid 
line in Fig. [1]). It has been speculated that collective neutrino flavor oscillations, during the 
post-bounce accretion phase, may affect neutrino luminosities and hence heating and cooling. 
Recently, it has been shown that matetr dominance suppresses collective flavor oscillations 
during the accretion phase [HI H] ■ It has been confirmed by several different groups based on 
different supernova models [8] [9] . 

At the onset of explosion, mass accretion vanishes and the central proto-neutron star con- 
tracts rapidly (see the magenta dash-dotted line in Fig. [TJ. It formed at core bounce and is 
hot and lepton-rich, in which terms it differs from the final supernova remnant neutron star. 
Between the expanding shock wave and the central proto-neutron star forms a region of low 
density and high entropy, where the surface of the proto-neutron star is subject to continued 
neutrino heating. There, a low-mass outflow develops known as neutrino-driven wind. The first 
sophisticated radiation-hydrodynamics study of the neutrino-driven wind was a milestone of 
research in the field [10]. It could explain the solar r-proccss abundances, due to the obtained 
strong wind with high entropies per baryon ~ 300 k# and generally neutron-rich conditions 
with Y e ~0. 35-0. 48. The neutrino-driven wind has also long been explored in static steady-state 
models [IT] Q21 Q2] [14] and dynamic studies [HI [TO] [16| [17] , as possible site for the nucleosyn- 
thesis of heavy elements [T5J [THl HO] . However, recent supernova simulations that include Boltz- 
mann neutrino transport cannot confirm the early results. They obtain generally proton-rich 
conditions and entropies per baryon on the order of 100 k^ [21] [22]. The main difference to 
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the early studies is related to the evolution of neutrino luminosities and mean energies. Within 
the current models, they reduce continuously during the proto-neutron star dclcptonization on 
timcscalcs on the order of 10 seconds after the onset of explosion. Furthermore, the v e and 
v e spectra become increasingly similar. Charged-current dominace reduces, because final state 
electrons become Pauli-blocked and nucleons become degenerate at the neutrinospheres, due 
to the increasing density. Instead, the spectra become dominated by neutral-current processes 
during the proto-neutron star dclcptonization. 

The relevance of collective neutrino flavor oscillations has long been investigated in various 
astrophysical applications. Although collective neutrino flavor oscillations are suppressed during 
the accretion phase, they may be relevant after the onset of explosion due to the continuously 
decreasing matter density in the presence of still high neutrino densities. Here, we explore the 
possibility of complete spectral flips of v e and v^/ T at a fixed flip energy and their impact to 
I'p-process nucleosynthesis for a selected trajectory from a supernova simulation of a massive 
iron-core progenitor. We find that it enhances the neutron production rate which in turn 
increases the production of heavy nuclei with A > 90. 

The manuscript is organized as follows. We will summarize main aspects of our core-collapse 
model in § 2. In § 3 we will illustrate standard neutrino-driven explosions of massive stars in 
spherical symmetry as well as the neutrino-driven wind phase after the onset of an explosion. 
§ 4 is devoted to the evolution of neutrino luminosities and spectra. Illustration of our simplified 
neutrino flavor flip analysis and the impact to nucleosynthesis will be discussed in § 5. We close 
with a summary in § 6. 

2 Core-collapse supernova model 

The simulations under investigation are based on general relativistic radiation hydrodynamics 
and three-flavor Boltzmann neutrino transport in spherical symmetry. For details, see the 
following references J23J EH E3 [26l EB ES ES] . Recent improvements of the adaptive mesh have 
been added in ref. [21]. It enables large stable timesteps and allows for long simulation times 
on the order of 10 seconds. The list of weak processes considered is given in Table [TJ including 
references. In addition, the implementation of the following weak process, v e + v e ^ v ^/r + ^V/r > 
has been discussed in ref [30], following ref. [31]. For the current study, weak magnetism 
corrections as well as iV-iV-recoil and ion-ion-corrclations have not been included in the weak 
processes. 



weak process 1 References 

1 u e + n±=,p + e~ [H] 

2 p e+ pt^ n + e + [32] 

3 v e + {A, Z - 1) ±5 (A, Z) + e" [32] 

4 u + N Uv' + N [32] 

5 v + (A, Z) v' + (A, Z) [32] 

6 v + e ± Uv' + e ± [31 H5] 

7 v + v ±5 e~ + e+ 01 [25] 
8 v + 9 + N + N^N + N [33] 



1 Notes: v = {v e , v e , v^/ T1 ^/ T }, N = {n,p} 
Table 1: Neutrino reactions considered, including references. 
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For matter in nuclear statistical equilibrium (NSE), the equation of state from ref. [34] was 
used. It is based on rclativistic mean field approach and the Thomas-Fermi approximation 
for heavy nuclei, with a simplified composition of neutrons, protons, a-particles and a single 
representative heavy nucleus with average atomic mass A and charge Z. Baryon contributions 
for matter in non-NSE are added using a slim nuclear reaction network for 20 nuclei (see [21] and 
references therein). It is used only for energy production. On top of the baryons, contributions 
from (e~, e + ) and photons as well as ion-ion-correlations for non-NSE are added [33] , 

Our core-collapse simulations are launched from the low-mass 8.8 M Q O-Ne-Mg-core [3"6"II3T] 
and from more massive iron-cores of 10.8, 15 and 18 M Q [35]. Their evolution during accretion 
and explosion, as well as the long-term evolution on timescales on the order of seconds after 
the onset of explosion, will be discussed in the next section. 




Radius [km] Radius [km] Radius [km] 

Figure 2: Radial profiles of selected quantities at the onset of explosion, for the 8.8 M Q O-Ne- 
Mg-core progenitor (data are taken form ref. [2"T]). 



3 Explosions and long-term evolution 

Neutrino-driven explosions of the low-mass 8.8 M Q O-Ne-Mg-core can be obtained even in 
spherically symmetric supernova models [31] [H] . The success of this model is related to the 
special structure of the progenitor. Only about 0.1 M of the 1.376 M Q core is composed of 
iron-group nuclei, at the onset of collapse. The outer layers are dominated by 20 Ne and 24 Mg 
as well as further out 12 C and 16 0. During collapse, the Ne and Mg layers are partly burned 
to iron-group elements and hence the enclosed mass inside the iron-core grows. Moreover, 
when the standing accretion shock reaches the interface between C-0 and He-layers, where 
the density drops over more than 10 orders of magnitude, it turns into a dynamic shock with 
positive velocities. It determines the onset of explosion, at about 35 ms post bounce, after 
which the shock expands continuously to larger radii. Fig. [2] illustrates the onset of explosion 
for this model (see ref. |21| and references therein). During the early explosion phase the P e and 
z^e-spectra are very similar. Note further, for this low-mass progenitor the shock wave expands 
basically into vacuum due to the extremely low density of the He-rich hydrogen envelop, where 
velocities on the order of the speed of light are reached. The competition between reactions 
(1) and (2) in Table [T] lead to even slightly neutron-rich conditions with Y e ~ 0.4681-0.4986 
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during the initial shock expansion after the onset of explosion between about 200-400 ms post 
bounce. The timescale for v e captures to turn material to the proton-rich side is not sufficient. 
In axially symmetric simulations, matter becomes even more neutron-rich early after the onset 
of explosion, developing mushroom-like pockets with Y e> m - m ~. 0.404 [33]. It may be a possible 
site for the weak r-process, producing elements with atomic mass between A >56-90, for which 
the production problem based on standard chemical evolution models has been discovered in 
ref. |44j . Recently, the question of an additional nucleosynthesis process required in order to 
explain the observed abundances of these elements has been addressed. It became known as 
light-clement primary process (LEPP) [35] and is an active subject of research. 




Radius [km] Radius [km] Radius [km] 



Figure 3: Radial profiles of selected quantities at the onset of explosion, for the 15 M Q iron-core 
progenitor. 

The situation is different for more massive stars, illustrated in Fig. |3j at the example of a 
15 M© iron-core progenitor. The extended high-density Si-S-layer surrounding the more massive 
iron-core, leads to a post-bounce accretion phase that can last for several 100 ms (depending 
on the progenitor model). The central proto- neutron star is much more compact at the onset of 
explosion. The luminosities of v e and 9 e are very similar during the post-bounce accretion phase 
as well as at the onset of explosion. Hence, due to the rest-mass difference between neutrons 
and protons, matter becomes proton-rich with Y e ~0.5-0.57. The magnitude of the differences 
between u e and v e luminosities and mean energies, and consequently Y e , is an active subject 
of research. It may change taking corrections from weak magnetism and improved weak rates 
into account. Note that the explosions for the iron-core progenitors under investigation are 
obtained applying enhanced heating and cooling rates (detailed balance is fulfilled), in order to 
trigger the explosions. For more details, see ref. [2JJ. Core-collapse simulations based on multi- 
dimensional models, that include sophisticated neutrino transport [30] [3TJ [3J] , are rec L u i re d for 
simulation times on the order of several seconds after the onset of explosion. 

After the onset of explosion a region of low density and high entropy develops between 
the expanding explosion shock and the central proto-neutron star. Moreover, at the surface of 
the proto-neutron star establishes net-heating, illustrated in Fig. |4(b)[ It leads to a low-mass 
outflow, known as the neutrino-driven wind. Compared to the very fast initial expansion of 
the 8.8 Mq model, the situation is different for more-massive iron-core progenitors. There, 
the neutrino-driven ejecta expand into the extended C-0 and He-layers with baryon densities 
between 10 1 g cm~ 3 to 10 3 g cm~ 3 , where also the shock expansion slows down. For illustration, 
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(a) Radial profiles of Y e and velocity. 
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(b) Heating rate during the neutrino-driven wind 



Figure 4: Electron fraction and velocity profiles at the early explosion phase at about 200 ms 



post bounce for the 8.8 M model in graph 4(a) and net energy-deposition rates at the onset 
of the neutrino-driven wind phase at about 1 second post bounce in graph [4(b) | (data are taken 
form ref. [2T|V 



see the radial profiles at a selected post-bounce time during the neutrino-driven wind phase in 
Fig. 5(a) at the example of the 10.8 Mq progenitor. Moreover, the neutrino-driven wind collides 
with the slower moving explosion shock. Note that in case of a super-sonic neutrino-driven 
wind, a reverse shock forms as shown in Fig. 5(b) This leads to an additional temperature 



and entropy increase (see Figs. 5(a) and 5(b) between 5000-6000 km). The impact of the 
reverse shock on possible nucleosynthesis has been investigated recently [IB]. At late times, 
neutrino heating reduces and the neutrino-driven wind turns back to sub-sonic velocities before 
it vanishes completely. 




Figure 5: Radial profiles of selected quantities during the neutrino-driven wind phase for the 
10.8 Mq progenitor model under investigation (data are taken form ref. [2Tj). 
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Figure 6: Post bounce evolu- 
tion of neutrino (dotted lines: 
v e , solid lines: D e , dashed 
lines: !/^ T , dash-dotted lines: 
v^/r) luminosities (top) and 
mean energies (bottom) for 
the 18 M Q progenitor model 
(data are taken form rcf. SI ). 
At the end of the accretion 
phase, indicated by the sharp 
jumps at about 350 ms post 
bounce, the luminosities of all 
flavors decrease continuously. 
The same holds for the mean 
energies. Furthermore, neutrino 
luminosities and spectra become 
increasingly similar for all fla- 
vors during the deleptonization 
phase. They converge for this 
model at about 20 seconds post 
bounce. It is related to the 
reducing dominance of charged- 
current reactions during the 
proto-neutron star deleptoniza- 
tion. Instead, the spectra are 
dominated by neutral-current 
reactions (neutrino-neutron 
scattering). 



4 Neutrino spectra evolution 

The evolution of neutrino luminosities and mean energies is shown in Fig. [5] at the example 
of the 18 M Q progenitor up to 22 seconds post bounce. The observables are sampled in the 
co- moving reference frame at a distance of 500 km, well outside the ncutrinospheres. 

The ^-luminosity, O(10 52 ) erg/s, rises slowly after the deleptonization burst has been 
launched at 20 ms post bounce. v e and v^/t are produced only after bounce. The v^/ T - 
luminosities rises until about 20 ms post bounce and the £/ e -luminosity rises continuously. After 
reaching their maximum, the i^/ r -luminosity decreases slowly continued during the accretion 
phase on timescales of 100 ms. Furthermore, the v e and v e luminosities are determined by mass 
accretion at the neutrinospheres. They rise slowly on a timescale of 100 ms and reach their 
maximum of several 10 52 erg/s at the onset of explosion at about 350 ms post bounce. After 
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that, mass accretion vanishes and the electron flavor luminosities decrease rapidly one order 
of magnitude within the first second after the onset of explosion. The i/^/ T -luminosity reduces 
accordingly and takes similar values as the electron flavor luminosities. The magnitude of the 
differences between the different flavors is an active subject of research. It depends sensitively 
on the weak processes considered and the dimensionality of the model. On a long timescale 
on the order of several 10 seconds, i.e. the proto-neutron star delcptonization, the neutrino 
luminosities of all flavors reduce below 10 50 erg/s. Furthermore, they become practically indis- 
tinguishable. In multi-dimensional models, and in the presence of aspherical explosions, mass 
accretion is still possible after the onset of an explosion. A possible enhancement of the neutrino 
fluxes remains to be shown for simulation times on the order of several seconds. 

The mean energies, shown at the bottom of Fig. [51 have a similar behavior as the neutrino 
fluxes. They rise during the early post bounce phase up to 12, 14 and 19 MeV for v e , v e and 
Vfj,/ T respectively, at the onset of explosion. After that, (E) u , decreases continuously. (E) Vc 
and {E)p c stay about constant until about 2 seconds post bounce, after which they decrease 
as well. The spectra of all flavors converge during the evolution after the onset of explosion, 
i.e. the difference between the mean energies of all flavors reduces continuously during the 
proto-neutron star delcptonization. For the 18 Mq progenitor model under investigation, the 
spectra have converged at about 20 seconds post bounce. It is related to the reduced domi- 
nance of charged-current reactions, due to final-state electron blocking and increasing nucleon 
degeneracy. Instead, the spectra are dominated by neutral-current reactions, in particular 
neutrino-neutron scattering. 

The small, and even reducing, difference between v e and v e luminosities and spectra has 
important consequences for the composition. It leads to generally proton-rich conditions with 
Y e ~0. 52-0. 56 for matter that becomes gravitationally unbound during the proto-neutron star 
delcptonization in the neutrino-driven wind. The results obtained for the low-mass O-Nc-Mg- 
core are in qualitative agreement with the results of the Garching group |22| . The magnitude 
of Y e obtained for the models under investigation depends sensitively on the equation of state 
and the weak processes used. 

5 Nucleosynthesis under proton-rich conditions 

Matter at the surface of the proto-neutron star is in NSE, due to the high temperatures and 
densities. During the expansion in the neutrino-driven wind, matter cools until reaching larger 
distance from the proto-neutron star surface, where nucleons recombine into heavy nuclei. This 
nucleosynthesis depends sensitively on the initial proton-to-baryon ratio. It is determined via 
the competition of reactions (1) and (2) from Table Q] in the dissociated regime at the proto- 
neutron star surface. It depends on the above discussed electron flavor luminosities and spectra. 
In the presence of similar v e and v e luminosities and mean energiefl matter becomes proton-rich 
due to the neutron-proton rest- mass difference |47j . The proton- rich conditions obtained lead 
to isospin symmetric nuclei, mainly 56 Ni, as well as 4 He and free protons. However, the further 
nucleosynthesis stops at, e.g., 64 Ge which has a long beta-decay half-life of ~ 64 s (known as 
waiting-point nucleus) and because 65 As has a low proton separation energy of ~ 90 keV. 
The situation changes including neutrino reactions, mostly v e because isospin symmetric 

For neutron-rich conditions, — 6v e < 4A, where e = (E 2 )/(E). (E) is the mean neutrino energy and 
(E 2 ) is the square value of the root-mean-square (rms) energy and the neutron-proton rest-mass difference 
A = 1.2935 MeV 
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Figure 7: Anti-neutrino spectra with (bottom) and without (top) spectral split in graph 7(a) 
and overproduction factors based on vp- process nucleosynthesis in graph [7(b) | (both figures are 
taken from ref. 1481 1. 



nuclei are inert to v e captures. v e can be captured at mainly protons as well as nuclei on 
timescales of seconds at distances of several 100 km when temperatures are as low as several 
10 9 K. With the consequently increased neutron density, waiting-point nuclei can be overcome 
via (ro,p)-reactions, and the mass flow can continue to heavier nuclei. This process is know 
as vp process [HJ [501 HO] • Which heaviest nuclei can be reached depends on the conditions 
obtained in the neutrino-driven wind. The full circles in the upper panel of Fig. |7(b)| illustrate 
the hnal abundances of a particular mass element from a simulation of a 15 M Q progenitor 
star [EQ [SI [50], labelled "1116 ms". The figure shows the ratio Af;/(A/ c jX ij0 ), where M l 
is the produced mass of isotope i with the corresponding solar mass fraction -X^.q and total 
mass ejected (taken form ref. [52]). The neutrino spectra required are taken from ref. [50j and 
approximated by an a- fit (for details, see [48] and references therein). 

Collective neutrino flavor oscillations have long been investigated in the context of core- 
collapse supernovae. Here, we explore possible effects of complete collective flavor flips between 



v e and ffi/T on the vp process, assuming normal mass hierarchy. Following ref. [531 . Fig. 7(a) 



illustrates the spectral flip taking place at energy of about 18 McV (bottom) in comparison to 
the unmodified spectra (top). Note that including nucleon recoil, especially for heavy- lepton 
neutrinos, will decrease high-energy spectral differences between v e and v^/ T an d hence the 
enhanced high-energy tail of the oscillated v e as illustrated in Fig. 7(a) will be reduced. 

For the vp process, v e spectra are required. In addition to the unmodified neutrino spectra, 
we include the oscillated spectra into the nucleosynthesis analysis and repeat the calculation. 
Relevant is the change of the neutron production rate due to the inclusion of the spectral 
flip. It increases the neutron production due to the enhanced high-energy tail of the flipped v e 
spectra. Furthermore, we parametrized the flip energy |48j . We found that the enhancement 
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of the neutron production rate is a robust result which is independent from the flip energy. It 
depends on the characteristics of the neutrino spectra obtained in particular simulations. Result 
of the nucleosynthesis outcome, including the complete spectral flip, is shown in Fig. |7(b)| (open 
circles). The neutron-production rate is enhanced by a factor of 1.4, according to the spectra 
taken from ref [52]. It significantly increases the production of heavy elements with A > 90. 
Furthermore, overproduction factors for nuclei with A = 64, 68 and 76 are reduced slightly. It 
has been discussed in more details in ref. |48j . 



6 Summary 

Core-collapse supernova simulations are investigated in spherical symmetry. We focused on 
the post-bounce accretion phase and the evolution after the onset of explosion, illustrating 
different conditions comparing the two intrinsically different core-collapse progenitors with an 
O-Ne-Mg-core and an iron-core. The explosion of the first one is a combination of neutrino 
heating and energy deposition from nuclear burning, on a short timescale of only few 10 ms 
post bounce. Furthermore, matter remains slightly neutron-rich for this progenitor only during 
the early explosion phase. On the other hand, the massive Si-layer surrounding the iron-core 
for more massive progenitors leads to an extended post-bounce accretion phase. It can last 
several 100 ms (depending on the progenitor model) during which the central proto-neutron 
star, which formed at core bounce, contracts continuously. The resulting similar electron flavor 
neutrino luminosities lead to generally proton-rich conditions at the onset of explosion. 

After an explosion has been launched, continued neutrino heating at the proto-neutron star 
surface leads to a low-mass outflow on a timescale on the order of seconds. It became known as 
neutrino-driven wind, for which we discussed and illustrated typical conditions. Furthermore, 
the neutrino luminosities and spectra of all flavors become increasingly similar during the long- 
term proto-neutron star deleptonization. It is a consequence of the decreasing importance of 
charged-current reactions, because at the neutrinospheres (a) electrons become Pauli-blocked, 
(b) the number of neutrinos reduces continuously and (c) nucleons become degenerate. Instead, 
the spectra are dominated by scattering at neutrons. The reducing difference between is e and v e 
leads to generally proton-rich conditions. It has important consequences for the nucleosynthesis 
of heavy elements and may allow for the vp process. 

Recently, it has been shown that neutrino oscillations are dominated by matter terms and 
hence collective flavor oscillations are suppressed during the accretion phase (6j El EJ [9j - However, 
they may be possible during the later evolution after the onset of an explosion, during which the 
matter density decreases continuously. In addition to the standard nucleosynthesis in proton- 
rich conditions, we explore the possibility of the vp process assuming a complete spectral flip 
between v e and v jl j r at a certain split energy. We used neutrino spectra from recent supernova 
simulations. It results in enhanced v e captures, due to the enhanced high-energy tail of v e , which 
in turn increases the number of neutrons present during i/p-process nucleosynthesis. It allows the 
matter flow to proceed to heavier nuclei and results in larger abundances of nuclei with A > 64. 
The current analysis has to be improved, using full energy- and angle-dependent neutrino 
oscillation techniques, in order to compute the spectral evolution consistently in massive star 
explosions. 
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